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摘 要 : 为 了 解决 蔬菜 运输 耗 时 长 、 成 本 高 、 保 鲜 时 间 短 ， 导 致 送 达到 客户 手 上 蔬菜 质量 降低 等 问题 ， 在 
考虑 了 车 辆 载重 和 时 间 窗 等 约束 条 件 下 ， 本 研究 提出 了 一 种 带 时间 窗 多 目标 蔬菜 配送 路 径 优化 的 遗传 -模拟 
退火 (Genetic Algorithm and Simulated Annealing, GA-SA) 算法 。 在 遗传 算法 (Genetic Algorithm, GA) 操 
作 过 程 中 引入 模拟 退火 (Simulated Annealing, SA) 算法 自 适 应 (Metropolis). 接受 准则 : 首先 将 原始 种 群 进 
行 遗 传 算法 的 选择 、 交 叉 、 变 异 等 操作 ， 形 成 新 一 代 路 径 种 群 ， 此 时 通过 引入 Metropolis 准则 ， 对 新 一 代 路 
径 种 群 分 布 情况 进行 修正 、 选 择 、 交 叉 、 变 异 ， 得 到 目标 路 径 种 群 ， 达 到 全 部 车 辆 配送 完 返 回 到 配送 中 心 
的 耗 时 最 少 、 成 本 最 低 、 和 车 辆 使 用 最 少 的 多 目标 ， 求 得 蔬菜 运输 的 最 优 路 径 。 设 计 以 保定 市 为 配送 中 心 以 
及 向 保定 市 下 辖 的 各 个 乡镇 为 配送 点 进行 蔬菜 运输 路 径 优 化 的 试验 ,结果 证 明 ， 与 传统 的 GA、SA 相 比 ， 
GA-SA 能 够 有 效 增 快 其 收敛 速度 ， 优 化 后 的 配送 路 线 总 成 本 分 别 降 低 了 约 23.7% 和 4%， 总 路 程 分 别 减少 了 
22.6% 和 3%， 耗 时 分 别 减少 了 26.2 和 2.6h， 车 辆 分 别 少 使 用 2 辆 和 1 辆 。 本 研究 可 为 冷 鲜 食 品 以 及 其 他 运输 
路 径 优化 研究 提供 参考 价值 。 
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1 引言 


近年 来 ， 随 着 人 们 对 新 鲜 蔬 菜 的 需求 日 益 增 
长 ， 中 国 设施 蔬菜 生产 发 展 迅 速 ， 已 成 为 世界 设 
施 蔬 菜 种 植 第 一 大 国 "。 但 蔬菜 是 易 腐烂 食物 ， 
对 运输 时 间 和 运输 环境 ”都 有 较 高 的 要 求 ， 需 
要 在 冷藏 状态 以 最 短 的 时 间 运 送 到 客户 手 上 ， 
此 运送 路 径 规 划 问 题 模 型 变 得 更 加 复杂 ， 是 典型 
的 非 确定 性 多 项 式 (Non-Deterministic Polynomi- 
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al, NP) 问题 3。 如 何在 一 定 的 约束 条 件 下 达 
到 成 本 最 低 、 耗 时 最 短 等 目的 ， 成 为 研究 者 重视 
的 研究 内 容 。Wang 和 Wasil ^ 对 近 五 十 年 来 关于 
研究 车 辆 路 径 问 题 进行 了 总 结 ， 从 精确 方法 和 局 
发 法 两 种 方法 解决 周期 性 车 辆 路 径 问 题 ( Vehi- 
cle Routing Problem, VRP) Jl [i] 3f p £ [o] 3t 
(Arc Routing Problem, ARP) 两 种 问题 进行 分 
析 , 但 是 只 考虑 了 常见 方法 。 周 洋 和 陆 轶 “ 对 A 
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蔬菜 运输 公司 物流 配送 路 径 进行 了 优化 研究 ， 采 
用 节约 里 程 法 对 其 进行 优化 ， 以 总 路 程 和 配送 成 
本 为 目标 ， 但 这 是 同一 种 货物 进行 配送 ， 没 有 考 
虑 物流 配送 中 的 各 种 相关 成 本 。 陶 涛 和 朱 家 明 TU 
提出 了 一 种 基于 Floyd 算 法 对 蔬菜 运输 路 径 进 行 
优化 ， 达 到 了 效率 最 大 化 ， 对 实际 生活 提供 了 一 
定 的 参考 价值 。 苏 贝 贝 “根据 应 急 物 资 配 送 过 
程 中 遇 到 的 公共 问题 以 及 自然 灾害 设计 了 路 径 规 
划 王 扰 管理 模型 ， 并 使 用 耦合 Dijstra 算 法 进行 求 
解 。Mendes 等 ”提出 了 一 种 在 线 降 维 方法 来 处 
理 有 需求 啊 应 运输 的 车 辆 路 径 问 题 。Lespay 和 
Suchan"" 针对 具有 时 间 窗 的 一 致 车 辆 路 径 问 题 
开发 了 启发 式 算法 。Molina 47 "^ 针对 具有 时 间 
窗 的 异 构 车 辆 路 径 问 题 ， 提 出 了 一 种 基于 蚁 群 系 
统 (Ant Colony System, ACS) 的 模 因 算法 。 赵 
邦 舌 ”改进 多 目标 蚁 群 算法 对 冷 链 物流 路 径 优 
化 进行 研究 ， 但 只 是 对 目标 进行 优化 ， 在 收敛 速 
度 上 尚未 改进 。 

上 述 研究 大 多 针对 单 目 标的 路 径 优化 ， 
并 且 在 收敛 速度 上 改进 较 少 "“， 选 择 的 是 少量 
客户 形成 的 小 数据 集 TS see f BEA In] T dc 
他 的 没有 时 间 窗 的 物流 运输 ， 也 不 同 于 软 时 间 
窗 ^U" 的 冷 链 运输 ， 不 同 的 蔬菜 有 不 同 的 保鲜 
度 ， 而 客户 又 需要 不 同 种 类 的 蔬菜 。 为 使 蔬菜 能 
在 保鲜 时 间 内 送 达 客户 ， 且 达到 配送 成 本 最 低 、 
耗 时 最 短 、 使 用 车 辆 最 少 的 目标 ， 本 研究 采用 三 
种 算法 进行 求解 。 考 虑 多 种 因素 ， 采 用 97 个 配 
送 地 点 和 1 个 配送 中 心 77, JE 98x98 AB E BSEC 
据 集 进行 仿真 ， 建 立 耗 时 最 少 、 成 本 最 低 、 车 辆 
使 用 最 少 的 多 目标 模型 。 并 针对 保定 市 以 及 周边 
乡镇 进行 蔬菜 配送 CU 进行 路 径 规 划 仿 真实 验 ， 
为 其 他 类 似 路 径 规划 77 提供 参考 价值 。 


2 蔬菜 运输 配送 模型 


车 辆 路 线 问 题 (Vehicle Routing Problem, 
VRP) 最 早 是 由 Dantzig 和 Ramser 于 1959 年 首次 
提出 所， 一 般 描述 为 一 定数 量 的 客户 ， 各 自 有 
不 同 数量 的 货物 需求 ， 配 送 中 心 向 客户 提供 货 


物 ， 由 一 个 车 队 负 责 分 送 货物 ， 组 织 适当 的 行车 
路 线 '”， 目 标 是 使 得 客户 的 需求 得 到 满足 ， 并 
能 在 一 定 的 约束 下 ， 达 到 路 程 最 短 、 成 本 最 小 、 
耗费 时 间 最 少 等 目的 。 

构建 蔬菜 配送 路 径 规划 模型 满足 以 下 约束 
条 件 : 

(1) 所 有 运输 蔬菜 车 辆 从 配送 中 心 出 发 ， 最 
后 回 到 配送 中 心 ; 

(2) 每 个 客户 被 一 辆 车 服务 ， 每 辆 车 走 一 条 
路 线 ; 

(3) 每 条 运输 线路 上 ， 客 户 需 求 量 不 能 超过 
车 辆 载重 总 和 ; 

(4) 每 条 配送 路 径 的 长 度 不 超过 配送 车 辆 一 
次 配送 的 最 大 行驶 距离 ; 

(5) 车 辆 不 能 有 重复 的 路 线 ; 

(6) 在 客户 需要 蔬菜 的 最 短 保 鲜 时 间 内 
送 达 。 

本 研究 的 蔬 沫 运输 路 径 规 划 问 题 可 以 描述 
为 : 多 辆 车 在 同一 个 配送 中 心 出 发 ， 分 别 向 多 个 
配送 地 点 运输 ， 每 辆 车 的 运载 能 力 一 定 ， 每 个 配 
送 点 的 需求 量 不 能 超过 每 辆 车 的 运载 能 力 ， 使 用 
同一 类 型 的 货车 ， 每 辆 车 的 保鲜 设置 和 封闭 性 都 
达到 相同 状态 ， 配 送 中 心 有 足 人 够 的 车 辆 被 调用 ， 
并 且 配 送 结束 之 后 返回 配送 中 心 ， 根 据 以 上 描述 
来 求 以 成 本 使 用 最 低 、 时 间 消 耗 最 短 、 和 车辆 使 用 
最 少 为 目标 的 最 佳 路 径 。 本 研究 的 案例 配送 地 点 
如 图 1 所 示 ， 其 中 黑色 圆 点 代表 配送 地 点 ， 三角 
代表 配送 中 心 。 

根据 蔬菜 配送 地 点 和 配送 中 心 创 建 数学 模 
型 ， 需 要 做 出 以 下 定义 : 共有 个 配送 点 ， 配 送 
车 辆 的 编号 为 c (c=1，2，…，7m)， 客 户 需求 量 
的 编号 为 @，(i=1，2，…，n)， 客 户 的 编码 为 i 
(i=1, 2, œ, n), EREN h, H m Až 
量 最 多 4 辆 。 因 此 建立 以 下 数学 模型 。 

目标 1: 成 本 消耗 最 少 ， 即 车 辆 运输 成 本 和 
车 辆 固定 使 用 成 本 的 总 和 最 少 。 

min Z, = SYM, Xn (1) 


i-0j-0k- 


其 中 ，Z, 表 示 消 耗 成 本 ,元 ; NN 为 配送 地 点 
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113.5 114.0 114.5 115.0 115.5 116.0 116.5 
经 度 /(*) 


图 1 配送 中 心 和 配送 地 点 位 置 点 图 


Fig. 1 Distribution center and distribution locations 


个 数 ， 个 ; K 为 车 辆 数 ， 辆 ; i、j 为 客户 编码 ; 
了 Yj 为 车 辆 从 客户 i 行驶 到 客户 j 的 行驶 过 程 ;M, 
为 客户 ;到 客户 7 所 消耗 的 成 本 ， 元 。 

目标 2: 蔬菜 运输 消耗 时 间 最 少 ， 即 每 辆 车 
行驶 的 总 路 程 除 以 蔬菜 运输 和 车辆 的 车 速 得 到 的 时 
间 Git oe E ED THEE FE INE TR] ILI EO, HL ZA 
7R, ho 

min, = Y v Yr, (2) 

其 中 ,为 车 辆 kx 从 客户 i 行驶 到 客户 j 的 耗 
时 ，h。 

目标 3: 使 用 车 辆 最 少 ， 用 乙 表 示 ， 辆 。 天 
表示 和 车辆 数目 ， 辆 。 

min Z, = K (3) 

为 了 达到 试验 目的 以 及 上 文 做 出 的 约束 条 
件 ， 其 约束 条 件 公式 如 下 : 


YYx<4 (4) 
c=1j=0 


AA <h (5) 
i=l 
Ys. =V (6) 
i-0 
3 = Yi (7) 
j-0 


N N K 
2.2/2. Te Tus (8) 


N N K 
be X T, (9) 
k-0 


i-0j-0k- 


其 中 ，4 需求 最 多 车 辆 数 ， 辆 ; c 为 车 辆 编 
码 ; /为 车 载 容 量 ，t; g 为 客户 需求 量 ，kg。 Tua 
为 客户 需求 蔬菜 的 最 长 保鲜 时 间 ，h; 7 蔬菜 开 
运输 时 间 ，h。 

在 上 述 数 学 模型 中 ， 公 式 (4) 表示 所 有 从 
配送 中 心 出 发 的 车 辆 数 在 给 客户 配 货 的 总 户 数 不 
超过 配送 中 心 自 有 的 车 辆 总 数 ; 公式 (5) 表示 
每 条 配送 路 线 的 总 客户 需求 量 不 大 于 车 辆 的 运载 
能 力 ; 公式 (6) 和 公式 (7) 表示 一 个 客户 只 能 
由 一 辆 车 进行 配送 ; AX (8) 表示 蔬菜 开始 配 
送 时 间 ; 公式 (9) 表示 满足 客户 需求 的 最 短 时 
间 内 送 达 。 


3 ”蔬菜 运输 配送 算法 设计 


3.1 遗传 算法 


遗传 算法 (Genetic Algorithm, GA) 属于 智 
能 仿生 学 算法 ， 思 想 源 于 生物 遗传 学 和 适 者 生存 
的 自然 规律 ， 是 按照 基因 遗传 学 原理 而 实现 的 一 
种 迭代 过 程 的 搜索 算法 所 。 本 研究 根据 配送 地 
点 确定 种 群 大 小 ， 之 后 进行 GA 的 选择 、 
变异 等 操作 。 

蔬菜 运输 配送 路 径 优 化 过 程 总 结 为 : 

Step 1: 编码 。 首 先 根据 蔬菜 运输 配送 路 径 
优化 的 特点 ， 根 据 编 码 的 生成 方式 ， 将 所 有 的 配 
送 地 点 自然 数 0、1、… 、N， 随 机 生成 全 排列 ， 
将 K-1 个 0 随机 插入 全 排列 中 ， 并 判断 每 条 路 径 
是 否 超载 ， 若 没有 超载 ， 为 有 效 染 色 体 ， 和 否则 为 
无 效 染色 体 ， 染 色 体 长 度 =N+K1。 

Step 2: 初始 化 群体 。 编 码 之 后 ， 随 机 产生 
1~N+K-1 个 个 体 ， 构 成 初始 种 群 。 

Step 3: 适应 度 函 数 。 判 断 需 求 点 7 的 需求 量 
&，kg， 是 否 小 于 或 等 于 和 车辆 的 载重 量 ， 寿 满足 
就 计算 相应 的 目标 函数 值 。 

Step 4: 选择 操作 。 遵 照 适应 度 越 高 ， 选 择 
概率 越 大 的 原则 ， 从 路 径 种 群 中 选择 两 个 个 体 作 
为 父 方 和 母 方 。 


n 
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Step 5: 交叉 操作 。 抽 取 父 母 双方 的 染色 体 ， 
进行 概率 为 P. (交叉 概率 ) 的 交叉 ,产生 有 效 的 
子 代 。 

Step 6: 变异 操作 。 为 保留 群体 中 个 体 的 多 
样 性 对 子 代 的 染色 体 进行 概率 为 P，( 变 异 概率 ) 


t 
变异 。 


Step 7: 重复 Step3，Step4，Step5，Step6 步 
又 ， 直 到 新 种 群 的 产生 。 


3.2 GA 改 进 


SA 算法 来 源 于 固体 退火 原理 ， 将 固体 加 温 
至 充分 高 ， 再 让 其 徐徐 冷却 。 加 温 时 ， 国 体内 部 
粒子 随 温 升 变 为 无 序 状 ， 内 能 增 大 ， 而 徐徐 冷却 
时 粒子 渐 趋 有 序 ， 在 每 个 温度 都 达到 平衡 态 ， 最 
后 在 常温 时 达到 基态 ， 内 能 减 为 最 小 。 

GA 易 过 早 收敛 是 效率 比 其 他 传统 方法 低 ， 
易 陷 入 局 部 最 优 解 ， 为 解决 这 一 问题 ,在 GA 的 
交叉 变异 等 操作 之 后 ， 引 入 Metropolis 准 则 ， 具 
体操 作为 : 将 原始 路 径 种 群 进 行 GA 的 选择 、 交 
又 、 变 异 等 操作 后 ， 形 成 新 一 代 路 径 种 群 ， 此 时 
引入 自 适 应 Metropolis 准则 ， 对 新 一 代 路 径 种 群 
分 情况 进行 修正 选择 交叉 变异 后 的 种 群 ， 得 到 目 
标 路 径 种 群 。 其 中 Metropolis 准则 常 表示 为 以 下 
模型 。 


1 E (X) «E (x T" 


E (x afe) -E EN 
aT 


exp E (X...) 2E (x T 


(10) 
其 中 ，P. 为 差 值 概率 ; aX Boltzmann W% 
引入 Metropolish 准 则 之 后 判断 是 否 保留 该 
个 体 的 过 程 如 下 述 (1) 和 (2) 所 写 。 
(1) 车 经 过 选择 交叉 变异 后 种 群 中 第 v 个 个 
体 的 能 量 值 ( 函 数值) 小 于 操作 前 第 4 个 个 体 的 
能 量 值 (也 数 值 )， 则 无 条 件 保留 该 个 体 ， 遍 历 
所 有 个 体 。 
(2) 若 新 个 体 v 的 能 量 值 CR AUT 
作 前 个 体能 量 值 〈 函 数值 ) ， 则 根据 两 者 状态 下 


的 差 值 概率 公式 进行 选择 优良 个 体 。 
^ 
p= exp[- 2 (1) 


公式 (11) 用 来 判断 是 否 选择 该 个 体 ， 即 车 
概率 为 Random [0, 1] 任意 数 ， 则 选择 遗传 操 
作 后 的 新 个 体 v; 反之 则 依然 保持 原始 状态 选择 
个 体 u， 此 为 Metropolis 准 则 。 

改进 后 的 GA 路 径 优 化 问题 的 算法 流程 如 下 
图 2 所 示 。 


目标 种 群 
图 2 遗传 -模拟 退火 算法 流程 图 


Fig. 2 Flow chart of genetic simulated annealing algorithm 
4 试验 设计 与 结果 分 析 


4.1 试验 数据 


试验 地 点 为 河北 省 保定 市 为 配送 中 心 ， 保 定 
市 下 辖 的 部 分 乡镇 为 配送 地 点 ,根据 各 个 配送 地 
点 的 实际 需求 、 保 定 市 盛产 的 蔬菜 品种 ”以 及 
每 种 蔬菜 的 易 腐 烂 程度 ， 把 蔬菜 分 为 以 下 三 大 
类 ， 如 表 1 所 示 。 

以 保定 市 为 配送 中 心 ， 向 保定 市 下 辖 的 部 分 
乡镇 进行 配送 。 试 验 包括 98 个 地 点 ， 其 中 包括 1 
个 配送 中 心 和 97 个 配送 地 点 ， 具 体 安排 如 附 表 
所 示 。 配 送 地 点 坐标 、 配 送 地 点 之 间 的 距离 等 数 
据 来 自 百 度 地 图 ， 采 用 人 工 采 集 的 实时 数据 和 根 
据 实际 情况 确定 。 数 据 的 处 理 与 编程 使 用 
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表 1 蔬菜 分 类 及 时 间 窗 


Table 1 Vegetable classification and time widow 


冷却 因子 为 0.99， 选 择 交 换 结构 的 概率 为 0.2， 
选择 道 结构 的 概率 为 0.5， 交 叉 概率 为 0.9， 变 异 
概率 为 0.05， 代 沟 概率 为 0.9。 根 据 实际 情况 确定 


蔬菜 种 类 ”适宜 温度 /*C 集装箱 温度 /*C 时 间 窗 /h ”蔬菜 分 类 
X 


0 8—12 0—4 一 类 


人 Me cn _ 类 车 辆 运输 成 本 为 0.5 元 km， 车 速 为 40 km/h， 配 送 
菌 类 蔬菜 0 8—12 0—6 三 类 车 的 车 载 容量 为 5 t。 

根茎 类 蔬菜 3—5 8—12 0—10 三 类 

瓜 果 类 蔬菜 78 8 0—10 三 类 4.2 试验 过 程 

2020aMATLAB 软件 完成 。 Step 1: 加 载 原 始 数据 ， 配 送 点 之 间 的 距离 


确定 算法 参数 ， 其 中 迭代 次 数 为 2000 次 ， ”单位 为 km， 如 表 2 所 示 。 
R2 蔬菜 配送 路 径 优化 试验 原始 数据 


Table 2 Original data of vegetable distribution route optimization test 


地 点 保定 市 定 州 市 涿州 市 e 腰 山 镇 WES seo ERLI 博 野 县 EH 
保定 市 0 
定 州 市 67.3 0 
涿州 市 99.3 158.7 0 

0 
TIL ERE 26.1 54.3 113.8 0 
浦 上 乡 33.6 50.6 112.4 8.8 0 
0 

北 果 

x$ 165.1 131.1 237,2 e 120.2 113.6 el 0 
博 野 县 49.1 53.7 141.1 er 100.1 95.4 er 156.9 0 

EH 60.3 66.5 142.9 e 91.8 87.6 er 173.6 14.8 0 


Step 2: 处 理 数 据 ， 给 每 一 个 配送 地 点 进 编码 ， 形 成 距离 矩阵 ， 单 位 为 km， 如 表 3 所 示 。 
表 3 蔬菜 配送 路 径 优化 试验 处 理 后 的 数据 


Table 3 Data after experimental processing of vegetable distribution route optimization 


地 点 编码 0 1 2 et 28 29 er 95 96 97 
0 0 67.3 99.3 e 26.1 33.6 e 165.1 49.1 60.3 
1 67.3 0 158.7 el 543 50.6 e 131.1 53.7 66.5 
2 99.3 158.7 0 e 113.8 112.4 EE 141.1 142.9 

0 
28 26.1 543 113.8 0 8.8 120.2 100.1 91.8 
29 33.6 50.6 112.4 8.8 0 113.6 95.4 87.6 
0 
95 165.1 131.1 237.2 el 120.2 113.6 EE 0 237.2 173.6 
96 49.1 53.7 141.1 el 100.1 95.4 e 156.9 0 14.8 
97 60.3 66.5 142.9 el 91.8 87.6 e 173.6 14.8 0 


Step 3: 随机 生成 路 径 初始 种 群 值 ， 车 辆 数 (Genetic Algorithm and Simulated Annealing，GA- 
目 为 22 辆 ， 和 车辆 行驶 总 距离 为 5798.6 km 。 SA) 三 种 算法 分 别 进行 优化 。 

Step 4: 使 用 GA、 模 拟 退 火 算法 (Simulat- Step 5: 得 到 配送 路 线 图 、 配 送 方案 路 线 图 

ed Annealing, SA) 以 及 遗传 -模拟 退火 算法 的 试验 图 和 总 路 程 、 总 耗 时 、 总 成 本 的 试验 
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结果 。 
4.3 结果 分 析 

以 蔬菜 运输 编码 的 一 次 运行 结果 为 例 进行 说 
明 。 对 算 例 进 行 了 2000 次 迭代 求解 后 得 到 了 图 3 ， 
即 本 研究 算法 运行 的 路 线 图 。 


39.6 


口 配送 地 点 
394r Liv 


113.5 114 114.5 115 115.5 116 116.5 
B/C ) 


图 3 蔬菜 运输 最 优 路 线 图 


Fig. 3 Optimal roadmap for vegetable transportation 


“0” 表 示 配 送 中 心 ， 数 字 “1 一 97” 表 示 配 送 
地 点 ， 车 辆 使 用 数目 为 20 辆 ， 车 辆 行驶 总 距离 为 
3564.79 km， 得 到 了 最 优 的 20 条 路 径 ， 如 表 4 
所 示 。 

4.3.1 算法 试验 结果 分 析 

fii FH Matlab 软件 进行 编码 ， 分 别 对 三 种 算法 
行 试验 经 过 编码 试验 ， 得 到 的 三 种 算法 最 优 
化 过 程 对 比 图 ， 如 图 4 所 示 。 

由 图 4 可 以 看 出 ，SA 收敛 过 程 波动 较 大 ， 
是 因为 在 SA 算法 过 程 中 冷却 速度 较 快 ， 收 敛 过 
程 无 法 稳定 ; GA 和 GA-SA 的 收敛 速度 趋势 基本 
相同 ,但 是 最 后 达到 的 最 优 解 不 同 ， 是 因为 传统 
的 GA 容易 过 早 收 敛 且 陷入 局 部 最 优 解 ; GA-SA 
能 够 跳出 局 部 最 优 解 ， 之 后 继续 训练 达到 最 优 
解 。GA-SA 在 迭代 次 数 200 次 左右 开始 收敛 ， 
最 后 达到 最 优 解 2000 元 ; GA 在 迭代 次 数 300 次 
时 开始 收敛 达 ， 最 后 到 最 优 解 2900 元 ; SATEXE 
代 次 数 1200 次 时 开始 收 僵 ,最 后 达到 最 优 解 
2000 元 。 可 以 看 出 ，GA-SA 的 收敛 速度 均 比 GA 
MSA, 更 快 的 找到 全 局 最 优 解 ， 且 GA-SA 求 


进 
进 


表 4 蔬菜 配送 路 线 
Table 4 Vegetable distribution routes 


路 线 编码 规划 路 径 地 点 编码 


1 0—228—223—531—525—553—226—50 


2 0—216—210—218—25—217—20 


3 0—566—535—234—^?0 


4 0 一 >27 一 >89 一 >90 一 >91 一 >95 一 >86 一 >94 一 >0 


5 0 一 >84 一 >87 一 >93 一 >92 一 >0 


6 0 一 >76 一 >41 一 >54 一 >21 一 >0 


7 0——70—53—5296—297—502—20 


8 0 一 >1 一 >81 一 >85 一 >49 一 >0 


9 0—>56 一 >2 一 >4 一 >22 一 >6 一 >7 一 >0 


10 0=>24=>30=332=3233 一 229 二 >0 


11 0—>61 一 >60 一 >59 一 >69 一 >0 


12 0—237—545—239—552—5240—551—250—^20 


13 0 一 >12 一 >14 一 >8 一 >9 一 >74 一 >0 


14 0—288—283—^282—20 


15 0—268—271—218—219—280—50 


16 0—220—219—-13—20 


17 0—572—264—211—263—261—251—258—^20 


18 Q—»73—2/5—265—215—511—20 


19 0 一 >30 一 >42 一 >38 一 >0 


20 0—>48 一 >47 一 >46 一 >44 一 >43 一 >55 一 >0 


得 解 比 GA 求 得 解 更 优 ， 虽然 最 后 SA 算法 能 
GA-SA 达 到 相同 的 解 ， 但 是 在 收敛 过 程 和 收敛 
速度 中 远 不 如 GA-SA, 证 明了 GA-SA 比 传统 的 
GA 和 SA 的 性 能 更 强 。 

图 5 是 GA、SA 和 GA-SA 收 敛 对 比 图 ，GA-SA 
算法 开始 收敛 的 平均 迭代 次 数 比 GA 算法 和 SA 
算法 分 别 少 118.4 和 117.5 次 。 即 GA-SA 比 GA 和 
SA 收敛 速度 快 、 效 果 好 。 

4.3.2 ”应 用 结果 分 析 

图 6 和 图 7 的 试验 结果 横 坐 标的 取 值 是 取 10 
次 实验 结果 的 平均 值 ， 分 别 是 GA、SA 和 GA- 
SA 在 相同 试验 条 件 下 车 辆 行驶 总 路 程 和 总 成 本 
的 比较 。 
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* SA — GASA == GA 


总 成 本 最 优 解 /元 


o 


0 500 1000 1500 2000 


迭代 次 数 /次 
图 4 GA、SA、GA-SA 最 优 解 进 化 过 程 示意 图 
Fig. 4 Schematic diagram of optimal solution evolution pro- 
cess of genetic algorithm, simulated annealing algorithm 


and genetic simulated annealing algorithm 


EH am sa EJ oe 


2000 


IGRA 


IRNAK 
图 5 GA、SA、GA-SA 收 全 对比 图 
Fig.5 Convergence comparison diagram of genetic algo- 
rithm, simulated annealing algorithm and genetic simulated 


annealing algorithm 


EB. M ES 


1 2 3 4 


5 6 
实验 次 数 /次 
图 6 车 辆 行驶 总 路 程 对 比 图 


Fig. 6 Comparison chart of total vehicle distances 


GA SA GA-SA 


1 2 3 4 5 6 ti 8 9 10 


实验 次 数 /次 


图 7 车 辆 行驶 运输 总 成 本 对 比 图 


Fig. 7 Comparison chart of total vehicle transportation costs 


由 图 6 可 以 看 出 ，GA-SA 配送 行驶 平均 总 路 
程 比 GA 和 SA 分 别 减少 22.6% (1044.64 km) 和 
3% (100.01 km)。 在 这 次 试验 结果 中 ，GA-SA 总 
路 程 最 大 值 为 3649.47 km， 最 小 值 为 3520.16 km, 
两 者 相差 129.31 km， 最 高 值 与 最 低 值 相差 较 少 ， 
证 明 试 验 过 程 优化 平稳 。 有 差 值 是 因为 GA 具有 
随机 性 ， 另 外 多 目标 优化 要 考虑 多 个 目标 之 间 的 
权衡 ， 很 可 能 是 为 了 使 车 辆 数 降低 ， 而 使 总 成 本 
的 值 变 大 。GA 总 路 程 最 高 值 达 到 4859.37 km, 
最 小 值 达 到 4417.9 km， 两 者 相差 441.7 km， 差 
值 较 大 ， 因 为 传统 的 GA 容易 过 早 收敛 。 从 而 陷 
人 局 部 最 优 解 无 法 跳出 。SA 算 法 总 路 程 最 高 值 
达到 3765.96 km， 最 低 值 达到 3624.29 km， 两 者 
相差 141.64 km。 总 路 程 无 论 是 最 大 数值 ， 最 小 
数值 ， 还 是 相差 数值 ， 本 研究 算法 的 试验 结果 均 
小 于 传统 的 GA 和 SA 算法 。 

蔬菜 配送 成 本 包括 人 工 成 本 、 机 器 损耗 1 
损耗 等 。 经 过 当地 调研 ， 纤 合 考虑 设 定 每 公里 消 
耗 成 本 为 0.5 元 ， 即 : 

总 成 本 = 总 里 程 X0.5 (12) 

可 知 总 成 本 和 总 路 程 是 正比 关系 ， 有 着 相同 
的 变化 的 趋势 。 由 图 7 可 以 得 出 GA-SA 算 法 的 配 
送 平均 总 成 本 比 GA 算 法 和 SA 算法 分 别 减少 
23.7% (553.1 元 ) 和 4% (62.26 元 )。 无 论 是 从 
最 大 数值 、 最 小 数值 ， 还 是 相差 数值 ，GA-SA 
总 成 本 的 试验 结果 也 都 小 于 传统 的 GA 和 SA 
算法 。 
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5 结论 


本 研究 对 叶 类 和 茄 果 类 蔬菜 、 菌 类 蔬菜 、 根 
茎 类 和 瓜 果 类 蔬菜 三 类 蔬菜 设置 三 个 不 同 的 硬 时 
间 窗 ， 一 类 蔬菜 的 时 间 窗 为 0~4h、 二 类 蔬 荣 的 
时 间 窗 为 0~6 h、 三 类 蔬菜 的 时 间 答 为 0~10h， 
对 蔬菜 运输 总 成 本 、 运 输 耗 费时 间 、 运 输 使 用 车 
辆 进行 多 目标 优化 ， 构 建 蔬菜 运输 配送 路 径 优 化 
模型 ， 并 对 GA 进行 改进 为 GA-SA， 求 得 蔬菜 运 
输 最 优 路 径 ， 得 到 结论 如 下 。 

(1) 本 研究 算法 有 更 快 的 收敛 速度 ， 比 GA 
和 SA 分 别 减少 118.4 和 117.5 次 ,更 短 的 迭代 
时 间 。 

(2) GA-SA 算 法 有 效 降 低 了 蔬菜 运输 总 路 
程 ， 同 时 在 蔬菜 保鲜 时 间 内 减少 了 运输 时 间 ， 降 
低 了 蔬菜 新 鲜 度 的 损耗 。 

(3) 本 研究 算法 有 效 降低 了 行驶 总 成 本 。 
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Abstract: There are higher requirements for the timeliness of vegetable transportation and distribution. In order to solve the prob- 
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lems of long transportation time, high total transportation cost and short preservation time of vegetables during transportation, 
considering the constraints such as vehicle load and time window, this study proposed a genetic simulated annealing algorithm 
(GA-SA) for multi-objective vegetable distribution path optimization with time windows. That was, the simulated annealing al- 
gorithm (SA) adaptive (Metropolis) acceptance criterion was introduced into the operation process of genetic algorithm (GA). 
The basic idea was: First, the original population was selected, crossed and mutated by genetic algorithm to form a new genera- 
tion of path population. At this time, by introducing metropolis acceptance criterion, and then, after modifying the sub situation 
of the new generation path population and selecting cross mutation, a new target path population was obtained. The improved al- 
gorithm retained the excellent individual, and the convergence speed, jumped out of the local optimal solution found based on 
genetic algorithm, and then found the global optimal solution. Then, the multi-objective of returning all vehicles to the distribu- 
tion center after distribution was the least time-consuming, the lowest cost and the least use of vehicles was achieved, and the 
optimal path of vegetable transportation was obtained. Taking Baoding city in Hebei province as the distribution center and 
some towns under the jurisdiction of Baoding city as the distribution points, the experiment of vegetable transportation path op- 
timization was designed. The experiments of genetic algorithm, simulated annealing algorithm and genetic simulated annealing 
algorithm were carried out, respectively. The comparative analysis was carried out from the aspects of convergence speed, total 
distance, total time, vehicles and total cost. The experimental results showed that, compared with the genetic algorithm and sim- 
ulated annealing algorithm, GA-SA could effectively accelerate its convergence speed. The total cost of the optimized distribu- 
tion route reduced by about 23.7% and 4% respectively, the total distance reduced by 22.6% and 3% respectively, the time con- 
sumption reduced by 26.2 and 2.6 hours respectively, and 2 and 1 vehicles were used less respectively. This study could also 
provide reference for the research of cold fresh food and other transportation path optimization. 

Key words: genetic algorithm; Metropolis guidelines; vehicle routing problem; vegetable transportation; simulated annealing al- 


gorithm; time consuming; cost; path optimization 
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附 表 :地 点 编码 表 


编码 地 点 编码 地 点 编码 地 点 “| 编码 了 点 编码 地 点 编码 地 点 编码 地 点 
0 保定 市 14 fume | 28 Æ | 42 长 古城 乡 | 56 泪水 县 70. EHS | 84 ”龙泉 关 镇 
1 定 州 市 15 留 村 乡 29 WES | 43 都 襄 乡 57 WE 71 RBS 85 PHH 
2 涿州 市 16 ” 正 村 乡 30 白云 乡 | 44 WERS | 58 清 殉 镇 72 MS 86 ”城南 庄 镇 
3 安国 市 17 BW 31 河口 乡 | 45 JUKS | 59 PEH 3 WW 87 KFX 
4 ”高 碑 店 市 | 18 户 木 乡 32 安阳 乡 | 46 罗 庄 乡 60 阳 城 镇 74 IE] BE 88 EWH 
5 易 县 19 ” 东 和 釜山 乡 | 33 台 鱼 乡 | 47 EK 61 魏 村 镇 75 ” 望 亭 乡 89 RIRS 
6 徐 水 县 | 20 KES | 34 唐 县 48 KFS 62 WEOE 76 WWE | 90 KRZ 
7 安 肃 镇 | 21 源 县 35 ”仁厚 镇 | 49 WI 63 KEH 77 BWE | 91 史家 寨 乡 
8 dE | 22 定 兴 县 36 王 京 镇 | 50. FZE | 64 大 庄 镇 78 安 新 县 | 92 Ves e 
9 大 因 镇 23 顺 平 县 37 高 昌 镇 | 51 羊角 乡 65 着 村 镇 79 雄 县 93 REOS 
10 — BIA | 24 WHA 38 北 罗 镇 | 52 石门 乡 66 白 团 乡 80 容 城 县 | 94 ”下 庄 乡 
11 清河 镇 | 25 KRS 39 HAH | 53 ”黄石 口 乡 | 67 ” 北 店 乡 81 曲阳 县 | 9%5 ” 北 果 元 乡 
12 ”高 林村 镇 | 26  ” 神 南 乡 40 ZWA | 54 倒 马 关 乡 | 68 AS 82 RFH | 96 WEE 
13 ”大 王 店 镇 | 27 ”高 于 铺 镇 | 41 川 里 镇 | 55 ”望都 县 69 FES 83 RFA | 97 Erg: 


